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Abstract 

Background: V\\\l-^ subtype B (HIV-IB) still dominates in resource-rich countries but increased migration contributes to 
changes in the global subtype distribution. Also, spread of non-B subtypes within such countries occurs. The trend of the 
subtype distribution from the beginning of the epidemic in the country has earlier not been reported in detail. Thus the 
primary objective of this study is to describe the temporal trend of the subtype distribution from the beginning of the HIV-1 
epidemic in Sweden over three decades. 

Methods: HIV-1 poi sequences from patients (n = 3967) diagnosed in Sweden 1983- 2012, corresponding to >40% of 
patients ever diagnosed, were re-subtyped using several automated bioinformatics tools. The temporal trends of subtypes 
and recombinants during three decades were described by a multinomial logistic regression model. 

Results: AW eleven group M HIV-1 subtypes and sub-subtypes (78%), 17 circulating recombinant forms (CRFs) (19%) and 32 
unique recombinants forms (URF) (3%) were identified. When all patients were analysed, there was an increase of newly 
diagnosed HIV-1 C (RR, 95%CI: 1.10, 1.06-1.14), recombinants (1.20, 1.17-1.24) and other pure subtypes (1.11, 1.07-1.16) over 
time compared to HIV-1 B. The same pattern was found when all patients infected in Sweden (n = 1 165) were analysed. Also, 
for MSM patients infected in Sweden (n = 921), recombinant forms and other pure subtypes increased. 

Significance: Sweden exhibits one of the most diverse subtype epidemics outside Africa. The increase of non-B subtypes is 
due to migration and to a spread among heterosexually infected patients and MSIVl within the country. This viral 
heterogeneity may become a hotspot for development of more diverse and complex recombinant forms if the epidemics 
converge. 



Citation: Neogi U, Haggblom A, Santacatterina M, Bratt G, Gisslen M, et al. (2014) Temporal Trends in the Swedish HIV-1 Epidemic: Increase in Non-B Subtypes 
and Recombinant Forms over Three Decades. PLoS ONE 9(6): e99390. dol:10.1371/journal.pone.0099390 

Editor: DImltrios Paraskevis, University of Athens, Medical School, Greece 

Received March 17, 2014; Accepted May 13, 2014; Published June 12, 2014 

Copyright: © 2014 Neogi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction In any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Anonymized sequences may be accessed via 
Genbank (accession numbers, JQ698667-JQ698874} [1 3]. In the present study, the sequences that were analyzed are representative of the entire country thereby, 
in principle, allow for the reconstruction of the transmission network. 

Funding: Funding provided by Swedish International Developing Agency, CHAIN FP7 EU, the Swedish Civil Contingencies Agency [SWE-2009-151], and the 
Swedish Research Council [521-2012-3476 and 2007-7092]. The funders had no role In study design, data collection and analysis, decision to publish, or 
preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: ujjwal.neogl@ki.se 

'9 These authors contributed equally to this work. 



Introduction 

Significant disparity is observed in the global distribution of 
human immunodeficiency virus type 1 (HIV-1) subtypes [1]. The 
greatest viral diversity is found in west and central Africa, while the 
epidemic is largely mono-phylogenetic in other regions, e.g. South 
Africa, Ethiopia, and India where HIV-1 subtype C (HIV-IC) 
dominates [2-5]. In resource rich areas, the epidemic was driven 
during the first decades by HIV-1 subtype B (HIV- IB) [6-8], but 



an increase of non-B subtypes in Western Europe has been 
reported [9,10]. However, cross-sectional data across 20 European 
countries show that HIV- IB is still predominant in all of the 
countries, ranging from 45% to 90% [7]. 

Representative information on the HIV- 1 subtype distribution is 
valuable for tracking the epidemiology and understanding the 
evolutionary trajectory of the virus. With the advance of 
automated subtyping tools, the accuracy of subtype prediction 
and the feasibility to analyse large datasets have increased. 
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However, subtyping of viral recombinant forms needs more 
stringent analysis and the use of atleast two of the bioinformatics 
tools available is recommended [11]. 

Sweden has a low prevalence (0.06%) and incidence of HIV- 1 
infection and the migration from high prevalence countries has 
therefore played an important role in the increased diversity of the 
subtypes since the 1990s [12]. Our previous study of patients 
diagnosed 1980-1993, showed that almost exclusively HIV-IB 
was circulating up to the end of the 1980s [12]. An increasing 
trend with non-B subtypes was thereafter observed with 30% of 
newly diagnosed cases in 1993 and we predicted that the HIV-1 
subtype distribution would become even more diverse [12]. A 
more recent study confirmed that HIV- IB has decreased over time 
in patients who were diagnosed between 2003 and 2010 in 
Sweden [13]. 

Although the HIV-1 epidemic, in Sweden thus shiftc-d from a 
mono-phylogenetic HIV- IB epidemic to a more diverse pattern 
over time [12,13], the detailed changes since the start of the 
epidemic up to now have not yet been described. We therefore re- 
subtyped all available sequences from HTV-l infected individuals 
diagnosed in Sweden between 1983 and 2012, and described the 
temporal trends the HIV- 1 subtypes and recombinants over three 
decades. 

Material and Methods 

Patients and dataset 

The data were derived from the Swedish InfCare HIV 
database, which at the 31" December 2012 included >99% of 
living adult Swedish residents with confirmed HIV infection 
(n = 6374), and the majority of all patients diagnosed 1983-2012 
(n = 9146), at 30 chnics from all regions of the country. The Public 
Health Agency of Sweden report 10332 individuals with HIV 
diagnosis during this time period (http://folkhalsomyndigheten. 
se/about-folkhalso myndigheten-the-public-health-agency-of-swe- 
den/) but in the official statistics the numbers are over reported 
due to duplicates. Of the 9146 individuals, 7068 (77%) subjects 
have been exposed to combined antiretroviral therapy (ART) with 
three drugs or more from 1996, and at 31'' of December 2012, 
5872 (92.1%) were active on ART. Since 1983, 747 patients have 
left the country, 61 have been reported missing and 1813 patients 
have died, as reported in the database. 

The database contains demographic data, clinical data, 
treatment history and the results of almost all routine genotypic 
resistance testing (GRT) ever performed, including the viral 
sequences, starting in 1992 [14]. In the national guidelines, routine 
GRT has been recommended for patients who fail ART and for 
newly diagnosed patients since 2001 [15]. In 587 out of the 2501 
(23%) patients who were diagnosed between 1981 and 1991, the 
viruses have been sequenced retrospectively. Of the remaining 
patients from the 1980s, >700 subjects have earlier been 
serotyped and/ or genotyped as belonging to mainly HIV- 1 B 
[12]. After excluding 33 HIV-2 patients, 253 patients with no 
registered date of HIV-1 diagnosis and 4925 patients with no 
sequence available in the database for re-subtyping, the present 
dataset included sequences from 3967 HTV-l infected individuals 
(Table 1). The year wise HIV-1 diagnosis among the different 
transmission risk group is presented in table SI. Among the 
patients without sequences diagnosed before 1992, 49.7% were 
men who have sex with men (MSM), 21.8% were intra-venous 
drug users (TVDU), 18.5% were heterosexuaUy infected and 10% 
were others/unknown. From 1992 up to 2012, the corresponding 
proportions were MSM: 20%; IVDU: 6.9%; heterosexuaUy 
infected: 58.2% and others/unknown: 14.9%. 



Sequences and subtyping 

For accurate identification of recombinant forms only partial /lo/ 
sequences that had a minimum length of 500 nucleotides (nt) were 
used. For patients with multiple sequences (1253/3967; 32%), the 
oldest sequence was used for subtyping and the most recent 
sequence was used to verify the subtype. 

HIV-1 subtyping was inferred initially by two methods: i) 
REGA subtyping tool version 3 (REGA v3), which uses an 
improved decision-tree algorithm and has greater sensitivity for 
identification of pure subtypes and recombinants [11]; and ii) 
COMET (http:/ /comet.retrovirology.lu/), which uses context- 
based modelling for expeditious typing of HIV-1 viruses. Both the 
tools showed best performances for pure subtypes and CRFs [1 1]. 
Additionally we used Recombinant identification program version 
3 (RIP 3.0) with 200 nt window size, which determines 
recombination in shorter serjuences (<700 nt) [http://www.hiv. 
lanl.gov/content/sequence/RIP/RIPexplain.html]; and HIV-1 
BLAST, available in the HIV Los Alamos Database [http:// 
www.hiv.lanl.gov/content/sequence/BASIC_BLAST/basic_blast. 
html], to identify the nearest HIV-1 sequences from difierent 
geographical proximity. Due to limitations such as differences in the 
length of the sequences, the presence of degenerate bases "N" inside 
the sequences (especially older sequences), thus detailed phyloge- 
netic analysis was not done. Further, the presence of drug resistance 
mutations can affect the phylogenetic inference but not affect 
subtyping and the length of the sequences could be adjusted in the 
alignment in the automated tools. However, the REGA v3 uses 
phyiogenetic methods to identify the subtype of a specific sequence. 

A sequence was designated as a particular subtype if two of the 
tiiree tools (REGA, COMET or RIP) determined it. The 
sequences were designated as pure subtype [A (Al and A2), B, 
C, D, F (Fl and F2), G, H, J, and K], circulating recombinant 
forms (CRFs) [CRFOl to CRF49 as of May 201,3], and unique 
recombinant forms (URF.s). The URFs were further characterised 
using bootscan analysis in[:orporated in the Rega v3 tool and 
validated in the SimPlot v3 software [16]. 

Statistical analysis 

Demography and clinical parameters were assessed by Kruskal- 
Wallis equality-of-populations rank t(;st for continuous variables 
and test for categorical variables. To assess trends of the HIV-l 
subtypes, a multinomial logistic regression model [17] was used for 
the 3967 patients with HIV diagnosis between 1983 and 2012. 
The outcome (the HIV-1 subtypes) was grouped into HIV- IB 
(reference group), HIV-IC, recombinants, and other pure 
subtypes. The year of HIV-1 diagnosis was the only covariate 
considered within the model. Square root of year of diagnosis was 
also added in the model to test for non-linearity, but was not 
significant. As sensitivity analysis and to get a more detailed 
description of which strains were primarily imported to Sweden 
and which strains actually have spread in Sweden, we performed 
analysc-s rc-strictcd to: 1) All j)aticnts infected in Sweden (n = 1551); 
2) Patients heterosexuaUy infected in Sweden (n = 307); 3) MSM 
infected in Sweden (n = 921); 4) All patients born and infected in 
Sweden (n=1165); 5) AU patients infected outside of Sweden 
(n = 2161); 6) AU patients born in Sweden and infected outside of 
Sweden (n = 455); 7) Patients born and infected outside of Sweden 
(n= 1699). Statistical analyses were performed using Stata version 
12.1 SE (StataCorp LP, USA) 

Ethical considerations 

The study was approved by Regional Ethics Committee 
Stockholm (2005/1167-31/3). The patient information was 
anonymised and de-identified prior to analysis. The authors 
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confirm that there are some restrictions on the data underlying the 
conclusions in the manuscript. Anonymized sequences may be 
accessed via Genbank (accession numbers, JQ698667^Q698874) 
[13]. The sequences that were analysed are representative of the 
entire country thereby, in principle, allow for the reconstruction of 
the transmission network. The entire dataset is available upon 
request from the steering committee of InfCare HIV. 

Results 

Patient dataset 

Patient's demographics are presented in Table 1. Self-reported 
country of infection was available for 91% (3600/3967) of the 
patients. Among them, 43% (1551/3600) of the transmission 
events had occurred in Sweden. The median year of diagnosis 
differed between the subtypes (HIV-IB: 2000; HIV-IC: 2005; 
other pure subtypes: 2005; recombinants: 2007) (p<0.001). 
Compared to HIV- IB, non-B subtypes had a lower median 
CD4-H T-cell count at diagnosis (B: 410/^1; C: 289/|xl; other pure: 
310/nl; recombinants: 330/|J,l) (p<0.001). Also, patients with non- 
B subtypes were more frequently late presenters than patients with 
HIV-IB (<200 cells/mm^: 27-34% versus 20%; <350 ceUs/mm^: 
53-62% versus 40%). 

Diversified presence of HIV-1 subtypes 

AU eleven pure subtypes and sub-subtypes within HIV-1 group 
M (3107/3967; 78%), 17 CRTs (757/3967; 19%), and 32 URFs 
(103/3967; 3%) were identified (Figure 1). HIV-IB dominated 
(47%) followed by HIV-IC (18%) and GRF01_AE (12%). HIV-IB 
dominated in MSM (91%) and IVDU (66%). Among FVDU, 27% 
were infected with CRF01_AE. A diverse pattern was seen in 
heterosexuaUy infected patients among whom HIV-IC dominated 
(3 !"'()). The city wise major pure subtypes HIV-IB and HIV-IC 
along with major CRFs 01_AE and 02_AG was presented in table 
S2. 

Subtype predictions using different tools 

The REGA v3 and COMET showed 96% concordance to 
determine the subtypes and recombinants forms. The majority of 
the discordance was observed between RIP 3.0 versus REGA v3 
and BLAST on subtyping of CRFs (mainly 01_AE and 02_AG). 
REGA v3 improved thus the identification of 02_AG compared to 
the RIP 3.0 tool, as reported earUer [1 1]. The highest concordance 
was observed between REGA v3 and BLAST (97%). For the 
patients for whom there was more than one sequence, a 98% 
(1207/1230) concordance was observed. The 2% discordance was 
due to read-length variations, mainly in the terminal region, which 
limited the possibility to identify recombinants. The median length 
of the sequences were 1058 nt (range: 551-1758 nt). 

Temporal trends of HIV-1 subtypes and recombinants 

There was a significant increase of newly diagnosed HIV-IC, 
recombinants, and other pure subtypes over time compared to 
HIV-IB (set as reference) (p<0.01) (Figure 2). The multinomicd 
regression model was adjusted for year of diagnosis tested for non- 
linearity however it was not significant. By extending the graph, 
the model predicts the exceed of recombinant forms by 2015. In 
the sensitivity analysis of year of HIV- 1 diagnosis, an increase of 
HIV-IC, recombinant forms and other pure subtypes, using HIV- 
IB as a reference, was seen for sdl patients (n= 1551) who were 
reported to be infected in Sweden. The relative risk ratio (RRR) 
for HIV-IC was 1.10 (95%CL 1.06-1.14), for recombinants 1.20 
(95%CL 1.17-1.24) and for other pure subtypes 1.11 (95%CL 
1.07-1.16) (Table 2). The same pattern was seen for patients born 



and infected in Sweden (n=1165) (HIV-IC: 1.06, 1.02-1.11; 
recombinants: 1.22, 1.17-1.27; other pure subtypes: 1.10, 1.0.5- 
1.17) and for patients heterosexuaUy infected in Sweden (n = 307) 
(HIV-IC: 1.12, 1.06-1.18; recombinants: 1.19, 1.1,3-1.27; other 
pure subtypes: 1.14, 1.08-1.20). For MSM patients infected in 
Sweden (n = 921), recombinant forms (1.27, 0.99-1.19) and other 
pure subtypes (1.1, 1.05-1.17) increased significantly compared to 
HIV-IB, but not HIV-IC (1.08, 0.99-1.19). 

In the second sensitivity analysis of the year of diagnosis 
restricted to the patients infected outside Sweden (Table 3), a 
significant increase of HIV-IC (1.07, 1.06-1.09), recombinant 
forms (1.15, 1.12-1.17) and other pure subtypes was seen (1.05; 
1.04—1.07) when all patients (n = 2161) were analysed as well as for 
patients (n= 1699) who were born and infected outside Sweden 
(HIV-IC: 1.08, 1.05-1.11; recombinant forms: 1.14, 1.11-1.17; 
other pure subt)'pes: 1.05, 1.0,3-1.07). For patients (n = 455) born 
in Sweden but infected outside, recombinants (1.16, 1.11-1.20) 
and other pure subtypes (1.08, 1.03-1.13) were significantiy 
increased, but not HIV-IC (1.03, 0.02-1.07). 

Mosaic pattern of URFs 

A total of 103 patients (3%) were infected by URFs. 
Recombinants between HIV-lAl and D (AID) were the most 
common (n = 21), foUowed by AlC (n=12), AIG (n=ll), BC 
(n= 10), BFl, BD, CD, and complex triple recombinants (n = 7, 
each) (Figure 3A). Further major URFs (AlC, AID, BC, AIG) 
which have HXB2 positions 2253 to 3260 and had minimum 
degenerate bases (^5) were used for bootscan analysis (n = 43) 
implemented in Rega v3 with 400 bp window size and 20 bp step 
size and SimPlot analysis with 300 bp window size and 20 bp step 
size. The breakpoint analysis showed 33 dilferent mosaic patterns 
of recombinants (Figure 3B), indicating a presence of diverse 
recombinants in the country. 

Movement of major HIV-1 subtypes and recombinant 

forms 

Based on the self-reported country of transmission along with 
the migration history in the country and probable geographic 
hnkages of the sequences using BLAST analysis, the predicted 
movements of the HIV-l strains were described. HIV-IC entered 
from three difiFerent regions: East Africa, Southern Africa and 
West-Central Africa (Burundi and Congo). CRF01_AE migrated 
from mainly Thailand and CRF02_AG from West-Central Africa. 
A fraction of 06_cpx migrated from Eastern Europe - Central 
Asia. HIV- 1 B was mainly introduced from Latin America, United 
States and Western Europe, although a small fraction was 
introduced from Thailand. 

Discussion 

The Swedish HIV-1 epidemic started in December 1979 in 
Stockholm with an outbreak of HIV-IB infection among MSM 
[18]. We describe now the changes of the diversity of the HIV-1 
subtypes during three decades, using viral sequences from nearly 
4000 patients, representing almost half of all diagnosed patients in 
Sweden ever. Our analysis shows that all eleven group M pure 
subtypes and sub-subtypes, 1 7 CRFs and several URFs have been 
introduced into the country and that there is a trend of increasing 
new diagnosis of HIV-IC and CRFs, as compared to HIV-IB. 

The initial spread of HIV-1 in Western Europe was driven by 
HIV-IB [19]. Although, HIV-IB was still predominant among tire 
MSM and the IVDUs in the present study, as also observed in a 
cross-sectional multinational European study [20], the prediction 
of a further diversification of HIV-1 subtypes in Sweden due 
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Total Heterosexual IVDU MSM 



■ Al ■A2 MB mC ID IFl ■F2 BG UH Mi 

■ K BOI AE ■02_AG ■ 03_AB n 06_cpx 107 60 ■ 09_cpx lO CD ll_cpx 12_BF 
13_cpx -18_cpx 20_BG 24_BG 25_cpx 33_01B 34_01B 37_cpx URF 

Figure 1. Distribution of HIV-1 subtypes in Sweden between 1983 and 2012. Final dataset included pol sequences (length >500 nt) from 
3967 HIV-1 infected patients. HIV-1 subtyping was determined by three methods (see method section). The sequences were designated as pure 
subtypes [A (Al and A2), B, C, D, F (F1 and F2), G, H, J, and K], circulating recombinant forms (01_AE, 02_AG, 03_AB, 06_cpx, 10_CD, 1 1_cpx, 12_BF, 
13_cpx, 18_cpx, 20_BG, 24_BG, 25_cpx, 33_01B, 34_01B, 37_cpx) and unique recombinant forms (URFs). 
doi:1 0.1 371 /journal.pone.0099390.g001 



migration was confirmed [12]. When virus sequences fi'om all 
patients were analysed, increasing trends of HIV-IC, recombi- 
nants (CRFs and URFs) and other pure subtypes were thus 
observed. This strong increase of viral diversity at the national 
level was to a large extent due to human mobility from high 



prevalence countries, both from those where such a viral 
heterogeneity pattern prevails and from countries where the 
epidemics are largely mono-phylogenetic, such as Ethiopia and 
Thailand. 




Figure 2. Proportions and predicted proportions of newly diagnosed HIV-1 subtype B (B), subtype C (C), recombinant forms (REC) 
and other subtypes (other) by year of diagnosis. Dotted lines represent actual prevalence in each year, while the smooth lines represent the 
predicted prevalence. Multinomial logistic regression model adjusted for year of diagnosis was used. 
doi:1 0.1 371/journal.pone.0099390.g002 
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In addition to the migration of already HIV-1 infected 
individuals, we found that patients who were born and infected 
in Sweden contributed to the changing HIV- 1 subtypes as well as 
those born abroad but reported to be infected in Sweden, as 
suggested by our sensitivity analysis. This was true also for MSM 
infected in Sweden where an increase of CRFs and other pure 
subtypes was described. Thus, the changing Swedish HIV-l 
subtype pattern is not only due to migration of persons who were 
infected at arrival, but also to a spread of non-subtype B strains 
among migrants after arrival to the country and MSM. It shall be 
noted that the country of infection was self-reported and data 
collected this way are likely to underestimate the proportion of 
migrants infected after arrival to a country [21]. It is therefore not 
unlikely that the infection with non-B subtypes among persons 
living in Sweden is even higher than reported in the present study. 
Our finding of an increase of HIV-IC infections among those 
reported to have become infected in Sweden is line with the 
reports that in other regions where HIV-lC has been introduced, 
its prevalence has increased gradually. E.g. in Brazil, a temporal 
increase in the prevalence of HIV-IC and related recombinants 
has been observed over time [22,23]. 

A contribution to the relative increase in non-HI V- 1 B infections 
among newly diagnosed patients in Sweden could be that 
undiagnosed HIV-1 infection is more common among migrants, 
as compared to MSM [24]. The time of diagnosis is shifted away 
from time of infection towards more recent years. Thus, the T-ceU 
counts at diagnosis were lower among non-B infected patients and 
late presentation was more common compared to HIV- IB 
infected subjects. Also, dififerences in the coverage rate and 
efficacy of ART may lead to a relative decrease in the number of 
persons infected by HIV- IB strains within the country since 
undetectable viremia due to ART is associated with minimal 
contagiousness [25]. Thus, although the general coverage rate in 
Sweden is presentiy high, >90%, the proportion of successfully 
treated MSM is even higher than among heterosexually infected 
individuals [24]. Considering the actual change in proportions 
over calendar years and the analysis predictions, the recombinant 
forms will probably dominate by the year 2015 among newly 
diagnosed patients in Sweden. But it is obvious that the 
development will still be highly dependent on the migration 
patterns for the coming years. 

There are some caveats that need to be mentioned. Half of the 
patients had not sequences available for subtyping, among whom 
MSM dominated among those diagnosed before 1992 and 
migrants among those diagnosed later. However, there was no 
dilference in gender and mode of transmission between those 
having sequences and those who did not have sequences. Also, 
most of these patients were diagnosed during the 1980s and we 
have earlier shown that the vast majority of patients during this 
time period were infected by HIV- IB [12]. At least 491 patients 
entered into the cohort when they were on therapy and due to 
non-detectable viral load, genotyping was not performed. How- 
ever, the majority of these patients were migrants from Africa. 
Altogether, the trend in increase of non-HIV- 1 B is unlikely to be 
over-estimated, rather the opposite. Another caveat is that the 
country of infection was assessed through self-reports, but this bias 
would also rather result in an underestimation of non-HIV- IB 
[21]. It shall however be noted that it is not uncommon the 
migrants visit their home country and that they might become 
infected there [21]. Thus, our study cannot precisely define the 
actual spread of non-B strains within Sweden. It is also likely that 
the prevalence of recombinant forms was underestimated since 
only the pol region was sequenced, as earlier described for the 
Indian HIV-1 epidemic when only one gene was analysed [3]. 
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Table 3. Patients infected outside of Sweden - sensitivity analysis of year of HIV-1 diagnosis and the trends In subtype pattern 
during 1983-2012, using subtype B (HIV-1 B) as reference. 



All patients infected outside 
Sweden 

(n = 2161) 



Patients born in Sweden and infected outside Patients born and infected outside 
Sweden Sweden 



(n = 455) 



(n = 1699) 



RRR (95% CI) 

HiV-1 B Ref. 
HiV-IC 

Year of diagnosis 1 .07 (1 .06-1 .09) 
Recombinants 

Year of diagnosis 1.15(1.12-1.17) 
Otfier 

Year of diagnosis 1 .05 (1 .04-1 .07) 



p-value RRR (95% CI) 
Ref. 

<0.001 1.03 (0.02-1.07) 

<0.001 1.16 (1.11-1.20) 

<0.001 1.08 (1.03-1.13) 



p-value 



<0.001 



RRR (95% CI) 
Ref. 

1.08 (1.05-1.1) 
1.14 (1.11-1.17) 
1.05 (1.03-1.07) 



p-value 



<0.001 



An univariate multinomial logistic regression model [1 7] was used. The outcome (the HIV-1 subtypes) was grouped into HIV-1 B (reference group), HIV-1 C, recombinants, 
and other pure subtypes. Relative risk ratio and 95% confidence intervals were determined. RRR represents relative risk ratio. 
doi:1 0.1 371 /journal.pone.0099390.t003 



Despite these caveats, we believe that the strength of the study is 
that it covers data from the beginning of HIV epidemic tUl date 
and includes almost all sequences performed in routine diagnostics 
in Sweden, corresponding to almost 4000 individuals. 



In conclusion, this is the first comprehensive study to describe 
the trends of the HIV- 1 subtype distribution in Sweden since the 
beginning of the pandemic. To best of our knowledge, the 
epidemic is one of the most diverse epidemics outside west-central 
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Figure 3. Types and mosaic pattern of HIV-1 unique recombinant forms identified in Sweden. A. Proportion of recombinant forms was 
identified based on REGAv3, COMET and RIP 3.0 followed by bootscan analysis incorporated in REGAvS and SimPlot analysis. B. Mosaic pattern of 
URFs identified by REGAvS and SimPlot. Major URFs (A1C, AID, BC, A1G) which have HXB2 positions 2253 to 3260 and had minimum 
degenerate bases (<5) were used for bootscan analysis {n = 43) implemented in Rega v3 with 400 bp window size and 20 bp step size and SimPlot 
analysis with 300 bp window size and 20 bp step size. The reference sequences were selected for SimPlot analysis is based on the recombinant 
pattern. Upper bar presents HXB2 co-ordinate. 
doi:1 0.1 371/journal.pone.0099390.g003 
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Afrk:a with presence of all pure HIV-1 M group subtype, 17 CRFs 
and 32 URFs. Countries such as Portugal, Belgium, Switzerland 
and more recently, UK and Germany are having similar trends 
with an increasing of non-B subtypes [20,26,27]. Presently, there 
seems to be no or limited transmission of HIV-1 between the 
major populations at risk but if sexual intermingling increases in 
the future, the viral heterogeneity may become a hotspot for 
diverse and complex URFs. Such a scenario is supported by the 
spreading of novel recombinants across multiple risk groups in the 
United Kingdom [28] and that several triple URFs already are 
circulating in Sweden. Apart from that HIV- 1 B was also found to 
be Continual molecular surveillance is recommended to under- 
stand the changing nature of the epidemic and the mobility of the 
HIV-1 strains. 

Supporting Information 

Table SI The Table describes the number of newly 
diagnosed patients whose virus was sequenced each 
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